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Abstract. Hints in the Fermi data for a 130 GeV gamma line from the galactic center have 
ignited interest in potential gamma line signatures of dark matter. Explanations of this 
line based on dark matter annihilation face a parametric tension since they often rely on 
large enhancements of loop-suppressed cross sections. In this paper, we pursue an alterna- 
tive possibility that dark matter gamma lines could arise from "semi-annihilation" among 
multiple dark sector states. The semi-annihilation reaction tl^itpj il^kl with a single final 
state photon is typically enhanced relative to ordinary annihilation '0^'^^ 77 into photon 
pairs. Semi-annihilation allows for a wide range of dark matter masses compared to the 
fixed mass value required by annihilation, opening the possibility to explain potential dark 
matter signatures at higher energies. The most striking prediction of semi-annihilation is the 
presence of multiple gamma lines, with as many as order N'^ lines possible for N dark sector 
states, allowing for dark sector spectroscopy. A smoking gun signature arises in the simplest 
case of degenerate dark matter, where a strong semi- annihilation line at 130 GeV would be 
accompanied by a weaker annihilation line at 173 GeV. As a proof of principle, we construct 
two explicit models of dark matter semi-annihilation, one based on non-Abelian vector dark 
matter and the other based on retrofitting Rayleigh dark matter. 
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1 Introduction 

While the gravitational evidence for dark matter (DM) is overwhelming [1-5], the precise 
nature of DM is still unknown. An attractive and well-studied possibility is for DM to be 
composed of long-lived neutral particles (WIMPs) [6-10]. In such scenarios, one typically 
expects feeble non-gravitational interactions between DM and standard model (SM) particles, 
especially if the relic density of DM is determined by thermal freeze out. Such interactions 
allow DM in the galactic halo to annihilate or decay to SM states, motivating a variety of 
DM indirect detection experiments. 

Recent hints of a 130 GeV gamma line in the Fermi public data [11-15] have highlighted 
the importance of galactic photons as a probe of DM. Since standard astrophysical processes 
do not yield monochromatic photons at these high energies,^ the presence of line emission 
from the galactic center is an unusually robust probe of DM [17-23]. The morphology of the 
gamma line feature is consistent with DM annihilation [24-28], in which case the measured 
photon energy is equal to the DM mass. Assuming a standard DM halo profile, the annihila- 
tion cross section is consistent with {av)^^^^^ 2^ 1.3 x 10~^^cm^s~^ [12], corresponding to 
around one-thirtieth of the expected thermal freeze out cross section. This tentative gamma 
line feature has motivated a number of new DM scenarios, based on both DM annihilations 
[29-44] and DM decays [45, 46], as well as constraints on such DM models from accompanying 
continuum photon emission [47-51]. While alternative interpretations of the Fermi line have 
been put forward [52-57], the DM hypothesis is still very appealing, and the DM community 
eagerly awaits an independent analysis by the Fermi collaboration. 

In this paper, we present an alternative interpretation of the 130 GeV gamma line 
as arising from "semi-annihilation" of DM [58] (see also Refs. [59-61]). Semi-annihilation 
appears in models with multiple DM species, and creates a gamma line through the process 

tpitpj ^ tpkl, (1.1) 
^See Ref. [16] for a discussion of potential astrophysical sources of gamma lines. 
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where '0^ corresponds to the i-th DM particle. Here, the monochromatic photon energy is 
given by the combination 



{nii + nij + mfc)(mj + mj — m^) 



(1.2) 



2{mi + nij) 



where rrii is the mass of ipi. In particular, if all DM species have the same mass, then 
a 130 GeV line corresponds to 173 GeV DM. More generally, semi-annihilation predicts a 
corresponding line for each ij k trio. 

Semi-annihilation possesses a number of features which make it interesting to study, 
even if the present Fermi hints ultimately have a non-DM explanation. 

• Parametrically larger cross section. For neutral DM, the annihilation process 
'^'0 77 is generated by loop diagrams with an a^y^ suppression factor,^ making it 
generically too small to account for the large {av)^^^^^ ^ 1.3 x 10"^'' cm^ s"-*^ cross 
section required to explain the 130 GeV feature. While semi-annihilation diagrams 
are also loop suppressed, '^i'^j il^kl is only proportional to a single factor of aEM, 
yielding a modest enhancement relative to the annihilation case. 

• Wide range of dark matter masses. In contrast to annihilation, gamma lines from 
semi-annihilation do not uniquely fix the DM mass. Therefore a signal at 130 GeV 
could arise from a wide range of DM masses. While not the focus of this work, this 
could allow a common DM explanation for both the 130 GeV gamma ray line and 
electron/positron excesses observed at higher energies [62-64]. Of course, one must 
account for the fact that at fixed DM density pdm = ^dm^dm, the indirect detection 
rate scales as n^y^ oc l/m^|^. 

• Accompanying annihilation signal. Generic models of semi-annihilation do predict 
an annihilation signal il^ii^i 77, albeit of weaker strength. In the case of degenerate 
DM states, the 130 GeV semi-annihilation line would be accompanied by a weaker 
173 GeV annihilation line, and detecting such a line would be a smoking gun for the 
semi-annihilation process. More generally, the El^^^ line should have E"^, eI^^ and E^ 
companions, offering the possibility of DM spectroscopy. 

• Generic absence of a 112 GeV feature. The ipip 77 annihilation diagram (with 

= 130 GeV) usually has a companion diagram '0'^ (with E^ ^ 112 GeV), and 

there are tentative indications of this second line in the data [13, 15, 24]. Indeed, if the 
112 GeV feature were not to persist, this would disfavor the annihilation interpretation 
since most annihilation scenarios yield two lines (see Ref. [24]). In the case of semi- 
annihilation, the companion diagram to il^iil^j ipkl is "^ii^j ipk^ which has no 
photon in the final state, and therefore no corresponding 112 GeV feature. Of course, 
there may be an ij k trio that "accidentally" yields a 112 GeV line, but this is not 
a robust prediction of the semi-annihilation framework. 

Given these unique features, it is worthwhile to consider various models and implications of 
semi- annihilat ion . 

The remainder of this paper is organized as follows. In Sec. 2, we review the semi- 
annihilation framework and explain the generic physics that leads to gamma lines. We 

^See however Ref. [44] for a tree-level explanation. 
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Figure 1: Two different spectra for the DM species adapted from Ref. [58]. In the case 
on the left, the three particles are mutually stable since their masses satisfy the triangle 
inequality. In the case on the right, the decay t/^k tl^itl^jc/)'' is kinematically allowed. 



present a new semi- annihilation model based on non-Abelian vector DM in Sec. 3, and show 
how to retrofit existing annihilation models to have enhanced semi-annihilation in Sec. 4. We 
conclude in Sec. 5. 

2 Gamma Lines and Semi- Annihilation 

The most well-studied models of WIMP DM involve a single relic particle stabilized by a Z2 
symmetry, as for example in supersymmetric models with i?-parity. However, any discrete or 
continuos global symmetry can be used for DM stabilization, opening the possibility for a far 
richer dark sector with multiple stable particles. Here, we focus on stabilization symmetries 
that allow for the semi-annihilation reaction, and consider the implications for indirect DM 
detection. We also briefly mention collider and direct detection expectations. 

2.1 Review of Semi- Annihilation 

For a generic dark sector with multiple stable species tpi, the semi- annihilation reaction 
reads [58] 

^pit/jj ^ tpk^l), (2.1) 

where are stable DM particles, and is either a SM field or a dark sector field which 
eventually decays to the SM. This reaction can take place if the DM is stabilized by a 
symmetry larger than just Z2.^ Such reactions are ubiquitous in multi-component models 
where the different species are stabilized by "baryon" and/or "flavor" symmetries, as in 
QCD-like theories. In semi- annihilation, the total DM number changes by only one unit, in 
contrast to ordinary annihilation where the total DM number changes by two units. 

The reaction in Eq. (2.1) could potentially make one of the DM species unstable, since 
the decay process tl^itl^jc/)* can be obtained by crossing symmetry. However, such decays 

may be kinematically forbidden. For example, in the absence of the weak interaction, the 
proton neutron n, and charged pion tt^ are mutually stable, and the process p7T~ nn^ 
with TT^ 77 would be an example of semi-annihilation. More generally, as long as the DM 
spectrum satisfies the triangle inequality rrik < rrii ^ ruj, as well as its crossed versions, the 
semi-annihilation reaction does not destabilize DM. This condition on the DM spectrum is 
depicted in Fig. 1. 

When semi- annihilation is present, DM dynamics in the early universe are far more 
varied than in the standard Z2 case [58-61, 65-67], especially in the presence of a matter- 
antimatter asymmetry [65]. In order to correctly compute the DM density, semi-annihilation 

"^The simplest example is a Z3 symmetry, where DM is composed of a single stable component with a 
semi-annihilation reaction ipip ipcj). 
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Figure 2: Schematic diagrams for sources of gamma ray lines from DM (semi-) annihilation. 
Loops of charged particles are depicted by the hatched boxes. We show: (a) ordinary annihi- 
lation, where the cross section scales as cr('0^'0^ 77) oc a^y^a^] (b) semi-annihilation, with 
cross section scaling as a{iljiiljj i/^k^f) ^ <^EM<^i<^j<^fc- 



must be included since it provides additional channels to dilute relic DM particles. Recently, 
the impact of semi-annihilations was implemented in micrOMEGAs software package for the 
cases of Zs and Z4 symmetries [66]. 

2.2 Gamma Lines for Indirect Direction 

Semi-annihilation not only affects DM production in the early universe but also leaves an 
impact today in DM indirect detection experiments. The differential energy flux of (f) particles 
from DM (semi-) annihilation in the Milky Way is of the form 

^r' ^^^^ (^^)^.-^ ^ - ^v) ■ (2-2) 

^ ij^k 

Here, the indices i and j run over the various different DM species with number density n^, 
and N'^J^^ (E^^^^) are the number (energy) of particles produced in each reaction with 
averaged cross section {av)ij^k- For non-relativistic DM in the initial state, the spectrum of 
(j) particles consists of several monochromatic lines with different intensities. The spectrum 
of observed SM particles in indirect detection experiments depends on the nature of (j). If the 
field (j) is unstable, we must to convolve the (j) spectrum in Eq. (2.2) with the decays of (j). 
Here, we focus on models where (semi-) annihilation directly produces gamma rays, namely 
</. = 7.4 

It is well-known that monochromatic photons can be produced in the DM annihilation 
process 77 with energy 

= m,. (2.3) 

A schematic diagram for this process is sketched in Fig. 2a. Given the very stringent limits 
on charged relic DM [68, 69], photons cannot be produced at tree-level through renormaliz- 
able DM-photon couplings. Instead, annihilation proceeds via loop effects involving charged 
particles (or corresponding higher-dimensional operators). Thus the cross section is expected 

^For the case of semi-annihilation resulting in a neutrino final state, = z^, see Ref. [67]. 
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to be proportional to which is why standard DM scenarios predict a feeble direct an- 

nihilation rate to photons. 

In contrast, gamma ray production from semi- annihilation scales with only one power 
of a^Mi as depicted in Fig. 2b. Moreover, dark sector couplings can easily exceed the electro- 
magnetic coupling while still remaining perturbative, so gamma ray production from semi- 
annihilations can be enhanced with respect to ordinary annihilations. 

As mentioned in the introduction, the semi-annihilation reaction tpit/jj tpkl leads to 
gamma lines at energies 

ij^k _ i^i + + mk){mi + ruj - ruk) . . 

~ 2K + m,) ' 

enriching the spectrum of gamma rays expected in DM indirect detection. In the simple case 
where all are equal, there are two lines at 

= (annihilation), E^^^^ — -m^ (semi-annihilation). (2.5) 

More general multi-component DM models with N species would yield N annihilation lines 
with energies given by Eq. (2.3), and as many as N{N — 1){N — 2)/2 semi-annihilation lines 
with energies from Eq. (2.4). Note that the resulting gamma ray spectrum from annihilations 
are directly related to the DM mass. On the other hand, the semi-annihilation spectrum is 
determined from the DM mass differences, and one could explain a 130 GeV line with all 
DM species in the multi-TeV mass range, or greater. Of course, the specific intensities of the 
lines will depend on the detailed couplings and symmetries of the DM scenario. 

2.3 Direct Detection and Collider Phenomenology 

For a given indirect detection signal, one must consult a specific model to determine the 
possible collider and direct detection implications. That said, we can make a few generic 
remarks about how semi-annihilation might modify DM phenomenology compared to ordi- 
nary annihilation. Usually, when one imposes the requirement of a particular annihilation 
cross section, either to explain the DM relic abundance or an indirect detection signal, it 
is possible to use the reversed Feynman diagram to predict collider production rates and 
signatures. Similarly, the same Feynman diagrams read on their side can be used to predict 
direct detection signals. 

In contrast, if the cosmology of DM is determined via semi- annihilation, this phe- 
nomenology is greatly changed. Collider production via the semi- annihilation diagrams would 
result in three DM particles in the final state, leading to a much smaller production cross sec- 
tion and modified kinematics. Similarly, direct detection via these diagrams would be quite 
suppressed, since elastic scattering diagrams only appear via DM loops. Of course, generic 
semi-annihilation scenarios often allow for standard pair annihilation, though from the para- 
metrics of Fig. 2, such annihilation diagrams are typically suppressed. For these reasons, 
one generically expects potential collider or direct detection signals from a semi-annihilation 
scenario to be suppressed compared to a corresponding annihilation scenario. 
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Figure 3: Gamma rays from semi-annihilation of vector DM through triangle diagrams. In 
general such processes are allowed, however if the vectors originate from a spontaneously 
broken gauge group and the charged fermions are vector-like under all symmetries, these 
diagrams are forbidden by anomaly cancelation. 

3 Semi- Annihilation of Vector Dark Matter 

We now construct an explicit class of models which give rise to gamma rays from semi- 
annihilation. Considering Fig. 2b, one could associate various spin-representations to '0^, 
and il^k' Particle assignments could include a pair of fermions and a boson, or three bosons. 
In supersymmetric scenarios, one can have particles of different spin transforming under a 
common DM stabilization symmetry. For simplicity, we restrict DM to consist of particles of 
the same spin (ruling out the participation of fermions), though we consider models involving 
fermions in Sec. 4. In this section, we focus on massive spin-1 vectors as our DM candidates. 
After introducing our vector DM model in Sees. 3.1 and 3.2, we discuss additional direct 
detection and relic abundance considerations in Sees. 3.3 and 3.4. 

3.1 Gamma Rays from Box Diagrams 

Vector DM can arise in a number of ways [59-61, 70-79], including as resonances of some 
strongly-coupled dark sector [60]. Here we focus on massive gauge bosons from a sponta- 
neously broken gauge group Gd^ which are stabilized by a remnant custodial symmetry [59]. 
For simplicity, we assume that the vector DM particles share a common mass My enforced 
by the custodial symmetry, though more general scenarios are possible. 

In order to generate gamma rays from semi-annihilation, we assume the existence of 
"messengers", namely matter which is charged both under electromagnet ism and Gd- We 
take the messengers to be fermions which carry vector-like charges under any symmetry. We 
also assume that the messenger masses Mm do not receive contributions from the higgsing of 
Gd- As a consequence, for a spontaneously broken gauge group, the requirement of anomaly 
cancellation forbids semi- annihilation through triangle diagrams such as Fig. 3. 

This leaves only the box diagrams of Fig. 4a as the origin of DM semi-annihilation. These 
diagrams have a kinematic structure analogous to the light-by-light scattering in QED [80- 
82]. The calculation of semi-annihilation through box diagrams is somewhat more involved 
than light-by-light scattering since three of the external bosons are massive. That said, 
general expressions for the QED box diagram were calculated for off-shell external photons 
by Karplus and Neuman in Refs. [83, 84], and these results can be employed for semi- 
annihilation. Extending this calculation to the non-Abelian case is straightforward. To 
lowest order in the messenger mass, the effective operator which is generated is a non-Abelian 
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(b) 

Figure 4: Production of gamma-ray lines from (a) semi-annihilation and (b) annihilation of 
dark gauge bosons for general dark gauge groups. Here, f are the generators and different 
gauge couplings are allowed for product gauge groups. Permutations of external legs are not 
shown. 



extension of the Euler-Heisenberg Lagrangian 
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(3.1) 



where the sum is over all possible 24 permutations of the four field strengths.^ This lowest- 
order effective coupling is sufficient for illustrative purposes, however we have calculated 
corrections at all orders in My /Mm by evaluating the Feynman parameter expressions found 
in Refs. [83, 84]. 

We will always assume that the messengers are more massive than the DM, i.e. My < 
Mmi forbidding direct DM annihilation into the charged fermions. Since the leading order 
term in the semi-annihilation amplitude scales as ^ (My/M^)^, we will require that the 
messengers are not too heavy. In the simplifying limit in which Mi — Mj = Mj^ = My and 
taking all dark gauge couplings to be equal, ^ the final result for 5-wave semi-annihilation 
{crv)s of species i and j into species k and a photon is 



1697 a^a^M 



4608007r 



M2 



MyV 
Mm) 



Mm \ 

My) 



(3.2) 



Here, F{Mm /My) is a form factor which parameterizes the deviation from the leading order 
result in an expansion in the vector to messenger mass ratio. Closed-form expressions in 
terms of Spence functions can be found by reducing the loop integral to the usual basis set of 
scalar loop integrals, however the final expression is extremely lengthly. We have calculated 
this form factor numerically and the result is plotted in Fig. 5. It is clear that for Mm ^ My, 
enhancements of 0{1) arise. 

The group-theoretic factor in Eq. (3.2) is 



t 



(3.3) 



^In the general case, the appropriate combinatorial factors should be included for identical fields, as in for 
photon-photon scattering. This is not a concern here, though, since the amplitude vanishes for any identical 
non-Abelian vectors. 

^The extension to product gauge groups is straightforward. 
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Figure 5: Enhancement of the full loop integral in comparison with the lowest-order terms 
found using the Euler-Heisenberg Lagrangian. For large messenger masses, the form factor 
approaches unity (dashed blue) and the Euler-Heisenberg Lagrangian is sufficient. As the 
messenger masses approach the vector masses, the higher-order corrections become impor- 
tant. The red dot corresponds to Mm/My — 200/173, which is the benchmark point used 
for the SU(3)rf model of Sec. 3.2. 

where the sum runs over the vector-like messengers transforming in the representation of 
Gd'^ Note that T^-^^ scales like Nf^ where Nf is the number of messengers, so {(7v)s scales 
like TV?. The total DM semi-annihilation cross section relevant to indirect detection is 



where fi is the fraction of DM made up by species i. With an unbroken custodial symmetry, 
fi = 1/Ny where Ny is the number of vector DM species. 

3.2 The Case of a Dark SU(3) 

To study the semi-annihilation process concretely, we now commit to a particular dark gauge 
group. While a U(l)^ gauge group could exhibit semi- annihilation, we know of no symmetry 
that would allow semi-annihilation without also allowing kinetic mixing between the dark 
U(l)(i and U(l)y, rendering the vector DM unstable. The case of SU(2)(i also does not work 
because SU(2)(^ has (pseudo-)real representations and thus T^-^^ = 0.^ The smallest SU(A^)(i 
group for which semi- annihilation can be achieved is SU(3)(i, and we will focus on this case. 

If the dark SU(3)(i is broken by a 3 x 3 bi-fundamental Higgs with vacuum expectation 
value proportional to the identity,^ then all of the SU(3)(i gauge bosons obtain the same 
mass, and are stable owing to the remaining custodial symmetry. In this case, all eight gauge 
bosons form equal components of the DM and each species makes up a fraction / = 1/8 

^Because the sum is over vector-like fermions, anomaly cancellation does not enforce T^-^^ — 0. Similarly, 
imposing a charge conjugation symmetry would not imply T^^^ = 0. 

One could consider the product dark group SU(2)cz x U(l)d, which has some elements with T^^^ / 0. 
^Such a spontaneous breaking pattern can be enforced with an appropriate flavor symmetry. 
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Figure 6: Contours in the Mm — <^d plane of constant ^{av)s{yV V^) for the SU(3)(i 
model with My — 173 GeV. These are given in units of 2.3 x 10~^^ cm^/sec, which is 
the approximate cross section required to explain the Fermi line after the DM density is 
adjusted to account for the increased DM mass. The form factor F{Mm /My) from Fig. 5 is 
included. For messenger masses Mm ^ My perturbative values of the dark gauge coupling 
are sufficient, and for higher masses one can see the expected oc scaling of the 

required coupling. 

of the DM density. Since the vectors all have the same mass, the vector masses must be 
My c:^ 173 GeV to produce a gamma line at 130 GeV. We will take Nf messengers in the 
fundamental of SU(3)(i carrying electric charge 1,^^ yielding the averaged group-theoretic 
factor 



Summing over all species in Eq. (3.4), the full 5- wave semi-annihilation cross section is 



where in the second line we have used the numerical value of the form factor for this choice 
of masses, and the cross section is divided by two since we only observe one photon per 
annihilation. This benchmark cross section is a factor of (173/130)^ larger than the nominal 
1.3 X 10~^^cm^/s value required for annihilating DM [12] to counteract the smaller number 
density for larger DM masses. The reference values were chosen to allow the desired semi- 
annihilation rate to explain the Fermi line using a quasi-perturbative value of < 47r, 
without requiring too much tuning between the messenger mass and the DM mass. The 
reference value Nf — 1 predicts the existence of three (= 3A^/) new charged fermions present 
just above the weak scale. Contours of {av)s{yV V^) in the Mm — oi^d plane are shown 
in Fig. 6, where the full form factor dependence on the messenger mass is included. 

^°We will discuss the possible SU(2)i, x U(1)y quantum numbers of the messengers in the next subsection. 
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Note that the rate for semi-annihilation is somewhat smaher than one would naively 
estimate from the diagram in Fig. 4a. The reason for this is twofold. First, the box diagram 
simply has a numerically small coefficient, beyond the estimate from naive loop counting. 
Second, there is a suppression coming from the large number of DM species (Ny = 8). Since 
there are only a small number of non- vanishing traces for semi- annihilation in Eq. (3.5), this 
leads to a fractional suppression of the effective cross section by l/Ny. Despite these sup- 
pressions, Eq. (3.6) show that semi- annihilation can explain the Fermi line for a perturbative 
value of a^. 

In addition to semi-annihilation, one has ordinary pair annihilation VV 77 shown 
in Fig. 4b. This yields a higher-energy gamma line at 2^ 173 GeV and two photons per 
annihilation, with a cross section given by 



where we have not included the form factor enhancement in this estimate. At present, a 
173 GeV line at this cross section is below current Fermi sensitivity, but would be a smoking 
gun for the semi- annihilation scenario if present. Finally, annihilation into 7Z is also possible, 
though the detailed rate depends on the electroweak charges of the messengers. Generically, 
the 7Z rate will be comparable to the 77 rate, yielding an additional gamma line at ^ 161 



It is clear that explaining the Fermi line via the semi-annihilation box diagram does 
require a large dark gauge coupling, with increasing values of needed for heavier messenger 
masses. On the other hand, Eq. (3.7) shows that an explanation of the Fermi line based on an 
ordinary annihilation box diagram would require an even larger gauge coupling (or a larger 
number of messengers). In this way, semi-annihilation processes can be enhanced compared 
to annihilation processes, making them a plausible DM explanation for the Fermi line. We 
will explore other ways that semi-annihilation can dominate over annihilation in Sec. 4. 

3.3 Messenger Direct Detection and Collider Bounds 

We have seen that SU(3)(i vector DM can yield a gamma line via semi- annihilation. Here, 
we consider possible constraints on such a scenario from direct detection and collider exper- 
iments. In the next subsection, we discuss relic abundance constraints. 

Without introducing additional couplings between the messengers and the visible sector, 
the messengers themselves are stable and thus contribute to the observed DM density. In 
addition, the messengers must consist of electroweak multiplets containing both charged 
and neutral fields, since if they only contained charged fields the model would face strong 
bounds on charged relics [68, 69]. For typical electroweak charges, radiative corrections split 
the degeneracy between the charged and neutral components, leading to a mass splitting 
M± — Mq ^ (9(100) MeV [85]. For this reason, charged relics are of no concern, since the 
charged fields will decay down to the neutral fields via virtual W boson emission. 

However, the neutral relics can still be a concern because of stringent direct detection 
constraints. For Dirac fermion messengers, the observed DM abundance is saturated for 
Mm ^ few TeV [85]. Since the abundance scales as Qm^^ then for Mm ^ 200 

GeV the messengers will comprise only an (9(1/100) fraction of the observed DM abundance. 




(3.7) 



GeV. 
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While this fraction is small enough not to affect the Fermi line, it is large enough to face 
direct detection bounds. Messengers with non-zero hypercharge experience spin-independent 
elastic scattering on nuclei mediated via Z-boson exchange, with a cross section well above 
current bounds even with a 1/100 dilution of the abundance. For this reason, we opt for 
SU(2)i, triplet messengers with zero hypercharge, since spin-independent messenger-nucleon 
scattering via Z-boson exchange is forbidden. In this case, a spin-independent cross section 
a 10~^^ cm^ is generated at the one-loop level [85], which is well below current bounds 
when combined with the low messenger abundance. 

Stable S\J{2)l triplet messengers also have an interesting collider phenomenology, sim- 
ilar to nearly pure winos in supersymmetric scenarios. The charged messengers can be pair 
produced weakly via the Drell-Yan process and then decay within a few centimeters to the 
neutral partner and a very soft pion. Such processes are extremely difficult to observe, 
however promising approaches have been demonstrated [86]. Detection of these messengers 
typically requires large integrated luminosities, and messengers with mass Mm ^ 200 GeV 
are consistent with current LHC bounds [87]. Depending on the lifetime of the charged 
messengers, potentially stronger bounds are provided by Ref. [88]. 

Of course, the main phenomenological feature of this work, namely gamma lines from 
DM semi- annihilation, does not rely on a specific model for messengers, and we have outlined 
only one possibility here. Similarly, if one allows for additional fields and couplings, then the 
messengers could be made unstable, removing any relic messengers and greatly changing 
the collider phenomenology. 

3.4 Relic Abundance 

A key question for any WIMP DM scenario is whether the DM relic abundance could be 
determined via thermal freeze out. As currently presented, the interactions in Sec. 3.2 are 
too feeble to yield the correct relic density. The loop-level (semi-) annihilation cross sections 
are too small to sufficiently dilute DM, and since Mm > My, tree-level annihilation to 
messengers is kinematically inaccessible at the desired freeze out temperature Tf 2^ My/20. 

There are a number of avenues to achieve the correct relic density. First, one can keep 
the current field content and lower the mass of the messengers such that My and Mm are 
nearly degenerate. In this scenario, DM freeze out tracks messenger freeze out until late 
times. This possibility has been discussed in Ref. [37] and can lead to the correct relic 
abundance via a form of assisted freeze out [89]. However, this requires a close coincidence 
between messenger and DM masses. 

A simple alternative is to employ the Higgs portal [90-94]. The fields which break 
the dark force are three fundamentals of SU(3)(i, and can be written as a 3 x 3 matrix of 
fields Hd- These fields can couple to the visible sector through the Higgs portal coupling 
C D X\Hd\'^\H\'^^ which respects the custodial symmetry. These dark Higgses obtain a vev 
{Hd) = Vdh/V^ and, after electroweak symmetry breaking, the diagonal component of 
mixes with the SM Higgs /iq, introducing couplings between the dark vectors and the SM 
Higgs. We can write these couplings as 

/:d5:1m^A-^;(i + ^ + ..Y, (3.8) 

^^Alternatively, we could mix the Dirac SU(2)i: doublets with a Majorana singlet to render the spin- 
independent process inelastic. 

^^As an example, messengers with non-trivial lepton number could decay to leptons and weak bosons. 
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Figure 7: Tree-level semi-annihilation of vector DM. This cross section scales with the 
mixing angle as {av)s{VV Vh) oc sin^ 9^ and can dominate freeze out even for small 
mixing angles. 



where 6h is the mixing angle between dark and SM Higgses, My = and the ellipsis 

denote additional dark Higgs states. Since the SM Higgs only mixes with the diagonal 
component of and the SU(3)(i generators are traceless, all cubic Vhh couplings vanish. 
This leaves the couplings VVh and VVhh, which are diagonal in SU(3)(i since Tr[t^t-^] = 6'^^ /2. 

The interactions in Eq. (3.8) allow vector DM to annihilate into pairs of SM fields and 
to semi-annihilate into a vector and a Higgs. Note that such processes are kinematically 
accessible because My — 173 GeV and rrih — 125 GeV. For simplicity, we will assume all 
other states (in particular the dark Higgs bosons) are heavier than the DM. 

We will calculate semi- annihilation first, since this will turn out to be the dominant 
interaction. The diagrams responsible for semi-annihilation are depicted in Fig. 7. Taking 
into account the relative fractional abundance of each vector, and summing over all species, 
the 5-wave DM semi-annihilation cross section is 

f^ (a 

dnal sin^ 9h 9M^ ) \ 

Inputting ruh = 125 GeV and My = 173 GeV yields 

1 VK) = 2.9 X (^)^ (^)^ (3.10) 

which is of the parametrically required size to explain the vector DM relic density [58].^^ 
For the expected Higgs decay channels, such thermal freeze out cross sections are below, 
but approaching, current indirect detection bounds on continuum photons [47, 48, 95],^^ and 
below current bounds on positrons, anti-protons [47], neutrinos [96], and radio emissions [97]. 

Since the dark gauge coupling is large, only a small mixing angle is required to maintain 
thermal equilibrium with the SM. The mixing angle must be less than sinOh % 0.0055, 
otherwise the vector DM abundance would be rapidly depleted. With this in mind, one can 
estimate the relative importance of other possible annihilation channels through their scaling 



^^Using the techniques of Ref. [58], one can find a semi-analytic solution for VLyh^ using the standard 
thermal freeze out approximation with the modified Boltzmann equation. Only one Boltzmann equation is 
needed, since the DM particles all have the same mass and carry the same DM number. The factor of 1/2 is 
important and can be understood as arising because each semi-annihilation only changes the total number of 
DM particles by one unit relative to annihilation. 

^^Each semi-annihilation typically results in a bb pair, each with an average energy of 76 GeV, and the 
relevant constraint can be roughly estimated from the constraints on 76 GeV DM annihilating to bb pairs, 
with the appropriate re-scaling due to different DM number densities. 
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Final State 



Scaling 



Vh 




hh 

/i* ^ WW, zz 
/l* tt 



adat sin^ 9h{p'^/M^) 



Table 1: Relative scaling of various DM annihilation final states. Here, /i* represents an 5- 
channel off-shell Higgs boson. With our benchmark values, the Vh semi-annihilation process 
dominates the DM relic computation. 

with the various couplings. In Table 1, we show the relative scaling for various annihilation 
final states. Since the mixing angle is small, annihilations to Higgs pairs will give at most a 
few percent correction. Similarly, annihilations to electroweak bosons through an 5-channel 
Higgs is also suppressed since ^ c^w- Annihilation to top quark pairs could lead to 
corrections as large as 10% at high temperatures since at ^ l/47r, however DM annihilations 
to top quarks is strongly phase space suppressed, leading to additional suppressions by a 
factor T/My ^ 1/20 at freeze out. Hence in this scenario, DM freeze out is dominated by 
semi-annihilation into a SM Higgs. 

As a final note, the required mixing angle is consistent with fits to Higgs observation 
data and also limits on additional Higgs-hke scalars [98]. Interestingly, if the dark Higgs 
states have mass below M^^ < 600 GeV then, due to mixing with the SM Higgs, they could 
show up at the LHC through their decays to WW and ZZ pairs. 

4 Retrofitting Existing Models: RayDM 

Going beyond the vector DM models already discussed, one could build models involving 
more general combinations of fermions, scalars, and vectors. The number of such possibilities 
is quite large, and we will not attempt an exhaustive classification. Here we focus on the 
interesting case of taking existing models that explain the 130 GeV line via annihilation, 
and retrofitting them to yield gamma lines from semi-annihilation. In particular, the large 
couplings that are typically required for annihilation can be reduced to more moderate values 
in semi-annihilation. In addition, the DM masses can take a broader range of values and still 
yield the 130 GeV feature. 

As an example, we consider a variant of Rayleigh DM (RayDM) [34, 42], where DM is 
either a Majorana or pseudo-Dirac fermion. These papers showed how the Fermi 130 GeV 
line could be explained by DM with mass ^130 GeV which annihilates through loops of 
messengers (electroweak charged scalars and fermions) to on-shell photons. For simplicity, we 
will consider a single Majorana fermion, in which case the transition dipole operator vanishes 
and all annihilation is achieved via the CP-odd Rayleigh operator. 

To retrofit RayDM models to incorporate semi-annihilation, we simply introduce an 
additional Z' which couples vectorially to the messengers. RayDM and retrofitted RayDM 
scenarios are depicted in Fig. 8. We choose this Z' to have a mass < Mz' < 2M^, such 
that the additional process XX ~^ 1^' is allowed but XX ~^ ^' ^' is kinematically forbidden. In 
general, one expects the Z' to kinetically mix with hypercharge and thus become unstable. 
While the process xX ~^ 1^' would then technically be an annihilation diagram since it 
changes DM number by two units, we will continue to call it a semi-annihilation diagram 
since it has semi-annihilation kinematics (i.e. and Mz' satisfy the triangle inequality in 
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'7 Xv 



'7 ^7 

(a) (b) 

Figure 8: Typical diagrams contributing to (a) RayDM and (b) a retrofitted semi- 
annihilation version of RayDM with an additional dark sector final state Z' . The loops 
involve both scalar and fermion messengers. 

Fig. 1). One could also envision semi- annihilation scenarios with a stable (or met a-st able) 
Z'-like state. 

The electroweak charges of the messengers are free parameters and so we must pa- 
rameterize this freedom. Attempting to remain as close as possible to the conventions of 
Refs. [34, 42], we adopt the convention that, after integrating out the messengers, param- 
eterizes the relative couplings of the DM to the hypercharge and SU(2)7, field strengths. We 
then define the additional parameter which gives the relative DM couplings to U(l)' and 
U(l)y. This yields the Rayleigh operators 



where the ellipsis denote additional couplings to B'B' which are unimportant for our pur- 
poses. As in Refs. [34, 42], the overall scale of the operator is 

' ^^'^ (4.2) 



A| 48M3^2 ' 

where A is the Yukawa coupling between DM and the messengers, and we have made the 
simplifying assumption that the scalar and fermion messenger masses are the same. 

Considering the usual annihilation scenario with a DM mass = 130 GeV, we find 
the annihilation cross section is approximately 

/^.A (^^ ^ ./.A - ^i^^(cQS CQS^ + sin sin^ Qwf ( M^V . . 



1.3 X 10-2^ cmVs(^)'( 



200 GeV \V M 



Mm J V 130 GeV 



X 



where we have chosen messengers with the electroweak quantum numbers of the Higgs, such 
that cos 9^ = 0.29. Hence quite a strongly-coupled, but still perturbative, theory is required 
unless the charged messenger mass is reduced close to the DM mass. 



^Also, for Majorana DM, CP-even couplings to XX would lead to p-wave suppressed annihilations. 
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Alternatively, with the addition of a Z', the Fermi line at 130 GeV could result from DM 
semi-annihilation with a larger DM mass. In this case, the semi-annihilation cross section is 



where M^/ = 2M^y^l — E^/M^ has been replaced such that — 130 GeV is ensured for 
any choice of M^}^ Again, this benchmark cross section is a factor of (173/130)^ larger than 
the nominal 1.3 x 10~^^ cm^/s value required for annihilating DM [12]. While a value = 10 
might seem on the large side, this is only due to the small size of the electroweak couplings, 
and = 10 simply corresponds to a U(l)' coupling 1. 

By using this retrofitting technique, semi- annihilation can extend more restricted anni- 
hilation models to allow for a wider range of DM masses and smaller couplings. It should be 
kept in mind that a line at 112 GeV from the 7Z final state is no longer a prediction of the 
model, though it can be easily accommodated if necessary by having multiple states. Also, 
we have not addressed the question of relic abundance, however the observed abundance can 
again be easily explained by extending the dark sector to include e.g. a dark singlet scalar 
mixed with the SM Higgs. 

5 Conclusions 

Regardless of whether the tentative 130 GeV Fermi gamma line persists, this intriguing 
feature is a reminder that monochromatic galactic photons offer a unique probe of DM, 
and it is therefore important to understand DM scenarios in which such gamma lines can 
be enhanced. In this paper, we have shown that semi-annihilation ipiipj tpkl can yield 
parametrically larger cross sections for gamma lines than annihilation tl^itl^i 77. This 
makes semi-annihilation a well-motivated target for future gamma line studies. Furthermore, 
unlike explanations of the the Fermi gamma line based on DM annihilation, semi-annihilation 
allows for a wide range of DM masses, presenting the opportunity for DM to impact indirect 
detection experiments at a variety of scales. 

We have highlighted two interesting semi-annihilation scenarios. First, we considered a 
vector DM scenario with a custodial symmetry which generates a 130 GeV line from semi- 
annihilation of 173 GeV DM. Finding a sub-dominant 173 GeV gamma line from annihilation 
to 77 would be a smoking gun for this scenario, as would proving the absence of a 7Z line 
at 112 GeV. In the context of an SU(3)(i dark gauge group, various group-theoretic factors 
conspire to require a larger than expected gauge coupling to explain the Fermi line, but there 
are indeed benchmark values of the messenger mass and dark gauge coupling where one still 
has perturbative control. 

Second, we have shown how existing annihilation models can be retrofitted into semi- 
annihilation models. Focusing on the case of RayDM, explaining the Fermi line with the 
annihilation diagram XX ~^ 77 ^he Rayleigh operator requires large (but still perturba- 
tive) Yukawa couplings to messengers. We have shown that including an additional massive 

^^This explains the linear scaling with M^, since the kinematic factor (1 — Mf//4M^)^ = (E^/M^)^. 




(4.4) 
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dark Z' allows the semi- annihilation process XX ~^ ^'l achieve the same 130 GeV feature 
for smaller couplings and a wider range of DM masses. 

As we eagerly await confirmation of the 130 GeV line, semi-annihilation is a reminder 
that multiple gamma lines offer the ability to perform DM spectroscopy. If the dark sector 
is sufficiently rich, then N species of DM can yield as many as 0{N^) lines through a 
combination of annihilation and semi-annihilation, with strengths determined by the various 
dark couplings and symmetries. Such gamma line measurements would then offer crucial 
insights into the nature of the dark sector. 
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